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Abstract. Using the density functional theory, we perform a full atomic relaxation of the bulk ferrite 
with 12.5%-concentration of monoatomic interstitial Cr periodically located at the edges of the bcc Fe^ 
cell. We show that structural relaxation in such artificially engineered alloys leads to significant atomic 
displacements and results in the formation of novel highly stable configurations with parallel chains of 
octahedrically arranged Fe. The enhanced magnetic polarization in the low-symmetry metallic state of 
this type of alloys can be externally controlled by additional inclusion of nonmagnetic impurities like 
nitrogen. We discuss possible applications of generated interstitial alloys in spintronic devices and propose 
to consider them as a basis of novel durable types of stainless steels. 



1 Introduction 

Last years demonstrate increased activities in the search 
for novel materials exhibiting controlled modification of 
electronic properties by inclusion or implantation of dif- 
ferent atoms or ionic groups. A prominent example of the 
implantation-altered systems is the stainless steel. In the 
steels, the implantation of chromium, molibdenium, ni- 
trogen and other chemical elements substantially changes 
the microstructure of subsurface layers and modify their 
corrosion resistance and hardness [1]. 

In the development of novel efficient multifunctional 
materials for technological applications in the long-term 
devices, the properties like hardness, corrosion, heat resis- 
tance and other types of mechanical and chemical dura- 
bility are of central interest [2,3]. It frequently appears in 
science and technology that well known materials doped 
by different chemical elements exhibit unexpected physi- 
cal properties not revealed previously. 

As an example of such a new unexpected behavior, in 
the present work we consider an alloy Fe-Cr. The alloys 
of Fe and Cr, doped by C, Ni and by other elements, are 
widely used as basic components for ferritic and marten- 
sitic steels. Substitutional alloys of Fe and Cr have at- 
tracted much attention of theory and experiment due to 
their rich magnetic properties characterized by local an- 
tiferromagnetism in the proximity of Cr atoms implanted 
into ferromagnetic iron [4,5,6,7]. Due to small differences 
between the atomic radii of iron and chromium, the mod- 
ification of the substitutional alloy properties is limited 



to the local magnetic transformation due to local changes 
in the electronic orbital occupancies, without significant 
structural modifications. In contrast to the substitutional 
structural configurations, the interstitial Fe-Cr alloys con- 
sidered in the present work contain Cr impurities which 
are located in the interstitial positions of the bcc lattice of 
Fcq;. In the recent theoretical studies of the Cr intersitials 
in Fe-Cr alloys, different types of interstitial configura- 
tions were analyzed. Among them, a pair configuration 
(111) dumbbell is considered as the most energetically 
favourable which requires about 4.2eV for its formation 
under irradiation [8,9]. 

In the present work, we consider a novel monoatomic 
interstitial configuration which contains single Cr atoms 
positioned in the centers of the edges of the bcc ferrite. In 
contrast to the substitutional alloys, the significant forces 
due to the interstitial atoms induce substantial structural 
optimization which enhances the volume due to modified 
lattice constants and leads to the relaxation of the atomic 
positions in the unit cell. We find that the relaxation of 
the initial bcc unit cell results in significant atomic dis- 
tortions and in the formation of atomic chain-like struc- 
tures. As appears in the density-functional-theory (DFT) 
calculations of the optimized structures, the energy gain 
achieved due to the structural relaxation of the considered 
interstitial alloy can approach 6.17 eV which makes this 
type of systems highly stable and durable. In the present 
work, we propose to consider these artificially generated 
alloys as candidates for novel types of stainless steels. 
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The fundamental difference between the industrial al- 
loys and the alloys studied in the present work is the 
ordered and periodic character of the latters. In the in- 
dustrial steels, the amorphic character of the systems is 
related to the random distribution of the impurities. The 
hardening of the steels proceeds through the surface treat- 
ment and is accompanied by formation of granular mi- 
crostructure with the spatially inhomogeneous impurity 
concentration and modified subsurface properties [11]. In 
the studies of the subsurface Cr-doped alloyed ferrite, we 
consider the supercells containing periodically located Cr 
atoms in the cubic lattice of Fcq,. The interstitial Cr in- 
duces significant atomic reconstruction with consequent 
break of initial cubic symmetry and stabilization of a new 
lower-symmetry state. The appearing structural transfor- 
mation has a character of a phase transition which occurs 
due to nanoscale tailoring of cubic Fe by interstitial inclu- 
sion of Cr atoms, the effect which can be experimentally 
verified by the means of modern methods like AFM spec- 
troscopy . 

Using the DFT-based structural optimization, we ob- 
tain the optimized atomic microstructure of a chain-like 
character where the chains of octahedrically arranged Fe 
atoms are formed along the (OOl)-axis. We find that the 
competing ferromagnetic and antiferromagnetic interac- 
tions lead to spatially inhomogeneous spin polarization. 
The magnetization of the structurally relaxed system is 
significantly enhanced as compared to the pure ferrite 
without Cr inclusions. The obtained enhancement makes 
the generated alloys perspective candidates for spin polar- 
izers in spintronic applications. In the generated chain-like 
structures, the relaxation is accompanied by the formation 
of spatial channels with extremely low carrier density. We 
suggest that these channels can be considered as paths for 
the low-barrier-migration of light impurities like H, N, Li 
or C. As an example of a light atom in the interstitial al- 
loy, we study of the migration paths of nonmagnetic nitro- 
gen and calculate the energy barriers along the migration 
paths. We obtain a strong influence of the nonmagnetic 
N on the alloy magnetization. Our findings show that the 
structural modifications due to possible nanoscale tuning 
of Cr impurities on the edges of bcc cubic cells of iron 
can play a central role in the control of their electronic 
properties. 

2 Structural relaxation of interstitial alloy 
Fe^-Cr 

The present studies of the electronic properties of the con- 
sidered interstitial Fe-Cr alloy are based on the DFT cal- 
culations of the electronic structure of the systems gener- 
ated by periodic translation of specially chosen supercells. 
The initial supercell shown in Fig. 1 contains the dou- 
bled 2x2 cubic bcc cell of ferrite (Fcq,) and a single Cr 
atom centered in one of the edges of the Fe^ cubic unit 
cell with the lattice constant a = 3.85 A. The obtained 
structure is described by a chemical formula FegCr and 
determines an interstitial Fe-Cr alloy with the Cr con- 
centration n = 0.125 which is typical for stainless steels. 



The presence of interstitial Cr leads to significant local 
forces acting on the neighbouring Fe atoms. To minimize 
the forces, the coordinates of all atoms have been relaxed. 
In the present studies, the optimization of the supercell 
has been performed by employing the DFT approach im- 
plemented withing the linearized augmented plane wave 
(LAPW) scheme in the full potential Wien2k code [13]. 
To study the role of the spin polarization in the structural 
relaxation, two different relaxation procedures have been 
employed. In the first procedure, the atomic optimal po- 
sitions are calculated in the local density approximation 
(LDA) on a 2 X 2 X 5 k-points grid. To explore the role 
of spin degrees of freedom in the relaxation, in the second 
procedure the local spin density approximation (LSDA) 
has been used in the optimization of the structure. The 
results of both methods of the structural relaxation are 
presented in Fig. 2. 




^ a 



Fig. 1. Schematic view of unrelaxed 2x2 Fe^ cell which 
contain 12.5% of edge-centered interstitial Cr. 



A central common feature which characterizes both 
(LDA- and LSDA-relaxed) structures is the clusterization 
of the sublattice of the iron atoms. In the LDA-optimized 
structure (Fig. 2(a)), the relaxation results in formation of 
a high-symmetry clusterized network. This network con- 
sists of the Fee-octahedra which form the square plaque- 
ttes in the (x, y) ((a, b)) plane with Cr atoms located in the 
center of each plaquette. The distance from the centered 
Cr to each nearest iron octaherda amounts 1.9 A. Despite 
the significant displacements of the iron atoms from their 
initial positions, the net electric polarization of the cell 
is zero due to high structural symmetry C4/m obtained 
after the relaxation. 

The formation energy of the relaxed FegCr-configuration 
can be expressed as 

Ef (LBA) = ^tot(Fe8Cr) - 8^tot(Fe) - ^tot(Cr), 
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where the last two terms identify the total energies of the 
bulk bcc Fea and Cr, respectively. To determine £'tot(Fe), 
we have calculated the total energy value of the bulk Fe^ 
in the ferromagnetic state. As the LSDA-calculation of the 
spin-polarized configurations of the bulk Cr are converged 
to the paramagnetic state, we consider the total energy 
£^tot(Cr) for the paramagnetic Cr. With these values, we 
find that Ef{LDA) = 4.82eV. To analyze the role of the 
relaxation, we have also calculated the energy Ef (umeV) 
of the formation of initial unrelaxed configuration which is 
equal to 5.02eV. As a consequence, the significant energy 
gain due to the structural relaxation 

Z\^(LDA) = ^/(unrel) - ^/(LDA) = 0.196ey, 

shows a central importance of the atomic displacements 
for the stability of the considered systems. 

The optimization procedure based on the LSDA ap- 
proach accounts for additional corrections due to spin po- 
larization and produces new ordered structural patterns 
presented in Fig. 2(b) and Fig. 2(c) for two different (un- 
relaxed a = b = 2. 86 A and relaxed a = b = 3 A) lattice 
constants. The volume-optimized structure (c) is signified 
by the 13%th increase of the unit cell volume due to the 
insertion of the interstitial Cr. The LSDA-optimized struc- 
tural pattern is characterized by the chains of atomic Fe- 
groups along the x(a)-direction, each group containing six 
Fe- atoms. The nearest chains are separated by a distance 
about 4 A and are connected to each other by the Fe-Cr 
bonds of the length about 2.4A for the structure (b) with 
a = 2. 86 A, and 2. 7 A for the structure (c) with the op- 
timized a = 3.0A. The local antiferromagnetic ordering 
in the vicinity of Cr is characterized by the magnetic mo- 
ments jiicr = —0.72 jiiB and /is = 2.4 jub and jhq = 1.25 /i^ 
of the neighbouring atoms Fe5 and Fe6, respectively. The 
magnetic moments of more distant iron atoms have the 
values around 2.5 /i^, which is close to results obtained 
for substitutional alloys and in pure Fe^ [8]. 

As compared to the tetragonal structure of the LDA- 
optimized system, the chain-like structure of the LSDA- 
relaxed supercell is characterized by substantially lower 
crystal symmetry and by the absence of the inversion 
center. In contrast to the LDA-based configurations, the 
formation energy of the LSDA-relaxed FegCr configura- 
tion Ef {LSD A) = Etot(Fe8Cr) - 8^tot(Fe) - ^tot(Cr) = 
— 1.15eV is negative which implies its high stability. We 
can also calculate the energy gain due to the structural 
relaxation by the LSDA approach 

AE{LSDA) = Efiuniel) - ^/(LSDA) = 6.17eV, 

which also demonstrates the high stability of the relaxed 
spin-polarized structure and a necessity to account for a 
spin polarization in the structural optimization of the sys- 
tems with strong magnetoelastic effect. 

2.1 Magnetic properties 

In the considered systems, we have also analyzed modifi- 
cation of the local magnetic properties due to the relax- 
ation of the interatomic distances. To see how the atomic 



displacements infiuence the spin polarization of the sur- 
rounding atoms, in Fig. 3 we present the dependences of 
the local moments of Cr and of two nearest neighboring 
Fe on the Cr displacement along the bond [Fe5-Cr-Fe6] 

A = [x{Ct) - x(Fe5)] - [x(Cr) - x(Fe5)]^ 

where [x{Ct) — x{Fe5)]^ is the optimized [Fe5-Cr]-bond 
length. The increase of A leads to the change of /icr from 
-0.7 to the value about -0.73 /j^b- In addition, the larger 
A implies the elongation of the [Fe5-Cr] bond and lead to 
the reduced /is = 2.39 jub due to the tendency for a sup- 
pression of antiferromagnetism in the vicinity of Fe5. The 
increase of A also produces an enhancement of /ie from 
1.25 jiiB to the values about 1.28 — 1.3 /i^, an opposite 
trend which occurs due to the shortening of the bond be- 
tween Cr and Fe6. 

In Fig. 3, the Z\-dependences of the atomic magnetic 
moments are highly asymmetric with respect to A. Conse- 
quently, the obtained magnetoelastic coupling produces an 
anisotropy of the magnetic moments and is accompanied 
by the loss of the inversion center due to the atomic dis- 
placements, the effect which can be observed in Fig. 2(b) 
and (c). In Fig. 2(c), the low-symmetry structure corre- 
sponds to the minimum of the total energy. As a con- 
clusion, the neglect of the magnetoelastic coupling in the 
electronic structure calculations does not allow to achieve 
a full optimization in this type of interstitial alloys. 

2.2 Electronic structure 

Fig. 4 shows the 3d spin-polarized electronic density con- 
tours of the LSDA-optimized structure in the (x, z) plane. 
One can see that the majority 3d spin- up states of Fe are 
highly occupied by the electrons whereas the electron con- 
centration of Cr spin up states is substantially lower. In 
contrast to this, the spin-down (minority) electrons are 
characterized by high electron occupation of Cr and lower 
electron density on Fe. In Fig. 4, the chain- like structures 
Fe-Cr in the 2:- direction are characterized by strong hy- 
bridization between the intra-chain 3d spin-down orbitals 
of Fe and Cr. The last feature leads to the spatial charge 
redistribution and to higher charge densities on the bonds 
between spin-down Cr and Fe. In the LSDA-optimized 
system, the structural optimization produces areas with 
low charge density in the y (6)-direction, where each area 
can be identified between the chains of Fe-octahedra. As 
can be seen in Fig. 4, these areas are almost free of the 
charge and can be considered as channels for the migration 
of light atoms like H, Li or N. Similarly to the contours 
in Fig. 4, the electron density of the majority Fe and Cr 
orbitals and on the bonds between Cr and Fe calculated 
for the LDA-optimized structure (Fig. 5) is substantially 
lower than the charge density on the spin-down contours, 
although the spatial charge distribution is more homoge- 
neous as compared to that in Fig. 4. 

For the LDA-relaxed structure, the density of states is 
characterized by strong suppression of the majority spin- 
up DOS at the Fermi level (Fig. 6(a)), whereas the minor- 
ity DOS at the Fermi level remains significant. Similar, 
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• (c) 



Fig. 2. Relaxed structure of Fe with 12.5% of Cr: (a) LDA 
calculations, (b) spin-polarized LSDA calculations in the 
structure with a = b = 2.86 A and (c) spin-polarized 
LSDA calculations in the structure with a = b = 3.0 A. 
The path 1 and path 2 identify possible pathes for diffusion 
through the channels formed due to atomic relaxation. 



although much stronger, suppression of majority DOS is 
typically observed in half-metallic systems where the elec- 
tric current is conducted by the electrons with the same 
direction of spin [14]. In contrast to the half- metallic- 
like features of the LDA-relaxed structure, the DOS of 
the LSDA-optimized system (Fig. 6(b)) demonstrates sub- 
stantial values at the Fermi level for both spin directions 
which implies an enhancement of the metallic state for the 
majority electrons. 

In transition metal oxides, the metallic state obtained 
in the LDA approach is usually strongly influenced by ad- 
ditional account for the local Coulomb corrections for the 




0.2 0.4 0.6 

change of distance [Cr-Fe5] 

Fig. 3. Local magnetic moments (in jj^b) of the atoms 
in Fe5-Cr-Fe6 triad versus the displacement zA=[Fe5-Cr]- 
[Fe5-Cr]0 of Cr along the (100) axis. Here [FeS-Cr]^ is the 
equilibrium distance between Fe6 and Cr. 



(i-electronic states [15,16,17,18]. In our work, the Coulomb 
corrections are incorporated within the SIC-variant of the 
LSDA+[/ approximation introduced in Ref. [15]. The re- 
sults are presented in Fig. 7 for two different values of 
U = 2 eV and U = 4.5 eV estimated and employed in 
Ref. [19,20,21] to account for the electron repulsion of 3d 
electrons of Fe and Cr. Fig. 7 shows the finite density 
of states at the Fermi {E = 0) level, although larger U 
leads to a significant suppression of the majority DOS at 
Ef which suggests a prevailing tendency towards a half- 
metallic behavior. 

In the LSDA-optimized structure, we find that the 
cell magnetic moment Mlsda = 3.84 fiB is larger then 
the magnetic moment Mlda = 2.88 jj^b in the LDA- 
optimized cell. Such an enhancement of the magnetic po- 
larization is connected with the substantial distortions 
ARi in the range 0.2-0.84 A and can be considered as a 
direct evidence of significant magnetoelastic effect. The lo- 
cal Coulomb corrections in the LSDA+/7-calculations re- 
sult in enhanced spin polarization. Specifically, we obtain 
Mlsda = 4.07 /ib for /7 = 2 eV, and Mlsda = 4.66 ^b 
for U = 4.5 eV. It is remarkable that in the substitutional 
alloy Fe-Cr with 12.5% of Cr, the LSDA approach gives 
the value 3.8 /hb for the cell magnetic moment which is 
slightly lower than the magnetic moment for the consid- 
ered LSDA-relaxed substitutional alloy. 

The obtained high spin polarization of the considered 
interstitial alloys Fe-Cr allows us to suggest these materi- 
als as possible candidates for spin polarizers in the spin- 
tronic devices. Possible technological applications of ar- 
tificially generated interstitial alloys would be related to 
the thin films produced for the needs of modern electronic 
industry. In such artificial systems, a central question is re- 
lated to the methods of implantation and positioning of Cr 
in centers of the edges of cubic bcc lattice of bulk Fe. With 
the current state of the art, such structural nanoscale ma- 
nipulation can be confined to the first subsurface layers of 
ferrite films by the using for instance the methods of opti- 
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Fig. 4. Contours of electron density maps in the (x, z)- 

plane {y/b = 0.2b^x and z given in A) obtained by integra- 
tion of electronic states in the energy window E between 
—3 eV below the Fermi level and the Fermi level. The 
results obtained by the structural optimization using the 
LSDA approximation. 



cal trapping by lasers [22] or by atomic force microscopy 
[11]. To estimate the stability of ultra-thin films of iron- 
chromium alloys which can be also considered as a basis 
for stainless steels, we have also extended our calculations 
to the nanoscale two-monolayers-thick iron films contain- 
ing one interstitial Cr atom per each 20-25 subsurface Fe 
atoms [12]. For such type of films, we have performed the 
calculations of the surface formation energy and of the 
electronic work functions, which were also compared to the 
corresponding quantities in the films of standard substitu- 
tional Fe-Cr alloys. In the interstitial Fe-Cr films, we find 
that the energy of the surface formation is about 1.69eV 
and the electronic work function amounts to leV, whereas 
for the substitutional Fe-Cr films we obtain 2.56eV for the 
surface energy and 0.57eV for the work function. These 
results allow us to expect high stability and durability of 
films generated on the basis of nanoscale-manipulated in- 
terstitial Fe-Cr alloys, as compared to the iron films with 
substitutional Cr impurities. 



spin up 




2 4 6 8 10 12 



X 

Fig. 5. Contours of electron density maps in the (x, z)- 
plane (y/b = 0.25, x and z given in A) calculated by in- 
tegration of electron states in the energy window E be- 
tween —3 eV below the Fermi level and the Fermi level. 
The LSDA results obtained in the initially LDA-relaxed 
structure. 

3 Migration paths of N in interstitial alloys 

Fec,-Cr 

A central question related to the stability of the consid- 
ered interstitial alloys is how various atomic impurities 
can modify the electronic properties. In the structurally 
relaxed alloys, the chains of atomic Fe-groups are sepa- 
rated by 4 A- wide atomic empty channels, which are ex- 
pected to contain pathways for light impurity atoms like 
H, N or Li. To explore a possibility of the migration of the 
impurities, we consider possible migration paths of a sin- 
gle nitrogen atom in the vicinity of the atomic Fe- chains 
in the interstitial FeQ,-Cr. 

In the studies of the migration paths of N impurities, 
we employed the nudged-elastic-band (NEB) method im- 
plemented in the Quantum-Espresso (QE) Package for the 
DFT calculations with the use of plane-wave basis sets 
and pseudopotentials [23,24]. In these calculations, for the 
atomic cores of Fe and Cr we employ the Perdew-Burke- 
Ernzerhof (PBE) norm-conserving pseudopotentials [25]. 
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Fig. 6. Total density of states (in eV~^) for structures 
optimized using (a) LDA approach and (b) spin-polarized 
LSDA approximation. The Fermi level corresponds to ^ = 
0. 



For each stage of the nitrogen transport, the NEB method 
involves a relaxation of the atomic positions and of the dis- 
tances between the different atoms in the supercell until 
the forces acting on the atoms reach their minima. In these 
calculations, we use the plane- wave cutoff 680 eV and the 
energy cutoff for charge and potential given by 1360 eV. 
In the NEB-approach, the relaxation of the atomic po- 
sitions along the nitrogen migration path is performed 
by the minimization of the total energy of each interme- 
diate configuration (image). These images correspond to 
different positions of N on the migration path and they 
are produced by the optimization of a specially generated 
object functional (action) with the consequent minimiza- 
tion of the spring forces perpendicular to the path. In 
our calculations, the convergence criteria for the norm of 
the force orthogonal to the path is achieved at the val- 
ues below 0.05eV/A. As the initial atomic configuration, 
the supercell Fe-Cr relaxed by the full-potential LSDA- 
approach [13] has been considered. 

Recent studies of the migration paths of single hydro- 
gen atoms by the pseudopotential NEB method demon- 
strate a good agreement of the obtained transport mecha- 
nisms and energy barriers with the experimental measure- 
ments [26]. Similarly to Ref. [26], in the present studies, 
each NEB-generated configuration has been modified by 
the introduction of the N atom and the obtained in this 
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Fig. 7. Total densities of states (in eV~^) for the LSDA- 
optimized structure calculated by the LSDA+/7 method 
with the local Coulomb corrections for the 3d-orbitals of 
Fe and Ct U = 2 eV (black curves) and 4 eV (blue curves). 
The Fermi level corresponds to E = 0. 



way extended supercell has been fully structurally opti- 
mized. 

To study the migration of N, we consider two different 
migration paths across the atomic Fe-chains in the Fe-Cr 
supercell. The first path (path (a)) describes the migration 
of N from the initial position inside the cell {z/c = 0.5) 
near the chain (1) across the channel to the chain (2) 
schematically presented in Fig. 8(a). In distinction to the 
path (a), the second path (b) reflects path2 in Fig. 2(c) 
and contains additional migration step of N from the su- 
percell boundary {z ^ 0) along the c-direction inside the 
supercell, with the further relocation through the channel 
to the atomic Fe-chain (2) indicated in Fig. 8(b). 

Fig. 9 shows the profiles of the total energy calculated 
along the N migration paths (a) and (b). The path (a) is 
characterized by the high energy barriers about 0.8 eV 
in the path coordinate range (0 < t/Rn < 0.4) and 
(0.8 < r/RN < 1) which corresponds to the migration of 
N within the two Fe-chains (1) and (2). Here Rn denotes 
the maximal length of the N path in the supercell which 
reaches about 1 nm for the path (a). The interchain mo- 
tion inside the channel is signified by a low energy barrier 
about 0.2 eV (0.4 < t/Rn < 0.7 in Fig. 9(a)). 

In contrast to this, the energy profile for the migra- 
tion path (b) (Fig. 9(b)) contains a plateau-like region at 
0.1 < r / Rn < 0.4. This miration step indicates the in- 
tracell replacement of N near the Fe-chain (1) along the 
[001]-direction demonstrated in Fig. 4, with a further relo- 
cation between the atomic chains across the atomic-empty 
channel with a low energy barrier about 0.2 eV. 

As a conclusion, we can note that the possible migra- 
tion paths of the light atoms in the considered intersti- 
tial Fe-Cr alloys contain a combination of the motion (i) 
within the atomic-empty channels and (ii) along the c- 
direction along to the Fe-contained atomic chains. 
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Fig. 8. Two different migration paths of nitrogen throu 
the channel of the optimized crystal cell of Fe-Cr. The t 
picture (a) represents the interchain migration of N insi 
the cell with the coordinate zj c near 0.5. The bottom p 
ture (b) corresponds to the migration of N from the c 
boundary [z = 0) along the z direction with the furtl 
migration between two neighbouring Fe-chains. The sy 
bols (1) and (2) denote the different atomic chains; (i) a 
(f) correspond to the initial and final positions of nitro^ 
in the migration paths. 



The question which arise due to the inclusion of N ii 
the magnetic Fe-Cr alloy is how the N impurities mc 
ify the magnetic properties of the system. In the work 
L Mazin[27], a comparison of the degrees of spin-polarize 
(DSP) calculated for Fe in the static limit through the d( 
sity of electronic states, via the current densities and 
the ballistic limit is presented. It is shown that all thi 
definitions of the DSP give very similar behavior for 
due to strong hybridization of the s'p and d states at t 
Fermi level. Thus we expect that in the considered Fe all^^^ 
with relatively low concentration of Cr it is sufficient to 
study the static spin polarization in order to capture the 
main properties of the alloy. Fig. 10 presents the change 
of the cell and atomic magnetic moments at the migra- 
tion of N along the path (a) and path (b). Although the 
nitrogen is initially nonmagnetic in the bulk, it becomes 
weakly magnetic inside the Cr-Fe alloy with a small mag- 
netic moment —0.04 \ib induced by the magnetism of the 
surrounding. It is noteworthy that the cell magnetic mo- 
ment is increased to 4 /i^ as N approaches Cr and the 
distance [N-Cr] becomes about 1.95 A. Such an enhance- 
ment of Mioi is explained by the strong atomic distortions 




N migration path coordinate (r/R^^) 



Fig. 9. Total energy profiles of the system along the mi- 
gration paths of N. Here denotes the energy of the 
system in the initial position of N and is the N coor- 
dinate in the final position of the path. 




N migration path coordinate r (A) 



Fig. 10. Cell (Mtot) and atomic magnetic moments (in 
\ib) along the migration paths of N. The red arrows iden- 
tify the maximal cell magnetic moments approached upon 
the minimization of the distance [N-Cr] along the N mi- 
gration paths. 



in the range between 0.04 A(Fe7) up to 0.2 A(Fe3) caused 
by the replacement of N and by the consequentmagnetoe- 
lastic effect. In Fig. 10, the increase of the distance from N 
to Cr suppresses the magnetic moment of N and decreases 
the cell magnetic polarization to the typical values about 
3.5 — 3.8 \iB obtained in LSDA-calculations for the ar- 
tificial Fe-Cr alloys. The obtained drastic change of the 
magnetic polarization clearly demonstrates a crucial im- 
portance of the location of nonmagnetic impurities like N 
for the electronic properties of alloy. As follows from our 
findings, a control of the location of N, for example by ex- 
ternal electric field, can lead to externally tuned changes 
of the magnetic polarization, a feature which is of central 
importance for possible spintronic devices based on the 
artificial Fe-Cr alloys. 



8 



N. Pavlenko, N. Shcherbovskikh, and Z.A. Duriagina: Title Suppressed Due to Excessive Length 



4 Conclusion 

We have performed the DFT studies of the bulk ferrite 
with 12.5%-concentration of monoatomic interstitial Cr 
periodically located at the edges of the bcc Fea cell. We 
have shown that the full atomic relaxation of the obtained 
interstitial Fe-Cr stabilizes a new chain-like low- symmetry 
structure. In this structure, the monoatomic Cr at the 
edges of ferrite bcc cells leads to the local atomic dis- 
tortions and results in the formation of parallel chains 
of Fee-ochahedra, which are connected by the interchain 
Fe-Cr bonds. The significant energy gain caused by such a 
structural relaxation approaches 6.17 eV which makes this 
type of interstitial alloy highly stable and energetically fa- 
vorable with the negative formation energy approaching 
— 1.15 eV. The novel electronic state of the system can be 
characterized as metallic, where the metallic properties is 
the result of strong Fe-Cr hybridization of the structurally 
relaxed alloy. In the investigations of the magnetic state 
of the generated relaxed structures, we have obtained a 
local antiferromagnetic order in the close proximity of Cr 
atoms, whereas the more distant Fe atoms are coupled fer- 
romagnetically. We also find that the nonmagnetic impu- 
rities like nitrogen can substantially modify the magnetic 
properties of the interstitial alloy which can be considered 
as an additional manifestation of the strong magnetoelas- 
tic eff"ect in this type of alloys. We propose to consider 
the generated interstitial alloys as perspective candidates 
for fabrication of novel highly durable stainless steels and 
for possible applications in spintronic and multifunctional 
devices. 
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